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ABSTRACT 


Principal  components  analysis  was  used  to  quantify 
subaerial  and  subaqueous  summer  to  winter  nearshore  morpho¬ 
logical  change  on  the  North  Carolina  coast  near  Duck.  The 
first  two  eigenvectors  explained  61.1%  of  the  total  vari¬ 
ance.  El  (35.4%  of  the  variance)  is  a  profile  rotation  and 
E2  (25.7%  of  the  variance)  is  a  bar/trough  function.  Two 
major  periods  of  change  were  found.  High  pressure  systems 
in  October  and  extratropical  cyclones  in  November  were  the 
weather  systems  implicated  in  these  changes.  These  storms 
caused  the  profiles  to  become  dissipative.  Steepness 
declined  and  two  bars  replaced  a  single  bar  profile  form.  A 
relationship  between  pre-existing  morphology  and  subsequent 
changes  in  the  equilibrium  bar/trough  morphology  was 
observed  during  the  storm  period.  Wave  characteristics 
(height,  period,  steepness,  direction,  and  breaker  type) 
could  not  be  associated  with  profile  changes.  The  simul¬ 
taneous  offshore  sediment  movement  across  three  well-defined 
zones  of  the  shorezone  suggests  that  shore-normal  edge  waves 
may  be  the  dominant  causative  agent  for  these  changes. 


INTRODUCTION 


Along  most  of  the  U.S.  coast  north  of  South  Carolina 
there  are  either  one  or  two  shore-parallel  bars.  Based  on 
the  analysis  of  color  infrared  aerial  imagery  along  the 
mid-Atlantic  coast  between  1970-1979,  Kochel  et  al.  (1983) 
determined  that  there  is  little  temporal  stability  in  the 
number  of  bars  present  from  Cape  Lookout  north.  Aerial 
photography  taken  over  the  last  decade  near  the  the  Army 
Corps  of  Engineers,  Field  Research  Facility,  Duck,  North 
Carolina,  showed  double  bars  normally  occurring  during 
winter,  while  only  a  single  bar  was  normally  present  during 
summer  months  (Table  1)  (Kochel  et  al.  1983,  1984).  The 
only  previous  mention  of  a  seasonal  transition  in  the  number 
of  bars  present  along  a  coast  comes  from  Shepard  (1950) .  Be 
observed  the  presence  of  a  single  bar  along  a  portion  of  the 
California  coast  in  winter  and  the  absence  of  a  bar  in 
summer.  Conventional  wisdom  (Evans,  1940)  and  theoretical 
studies  (Lau  and  Travis,  1973;  Carter,  Liu  and  Mei,  1973) 
indicate  that  the  number  of  bars  present  is  a  function  of 
the  overall  slope  of  the  nearshore  bottom  -  the  more  gradual 
the  slope,  the  greater  number  of  bars. 

Empirical  studies  using  multivariate  statistical  anal¬ 
yses  of  inshore  profiles  clearly  indicate  a  statistical 
independence  between  the  slope  of  the  profile,  the  curvilin¬ 
ear  form  of  the  profile,  and  the  number  of  bars  in  the 
profile  (Hayden  et  al.,  1975,  Resio  et  al.,  1974,  Dolan  et 
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al.,  1982).  It  is  unlikely  that  slope  variations  at  a  site 
are  responsible  for  the  temporal  instability  of  bar  numbers 
along  the  east  coast.  Our  understanding  of  the  nature  and 
causes  of  changes  in  nearshore  bathymetry  is  based  largely 
on  laboratory  and  theoretical  studies.  Spatially  and 
temporally  extensive,  systematic,  and  ongoing  beach  obser¬ 
vational  programs  have  only  recently  begun  at  Duck,  North 
Carolina. 

The  summer  to  winter  morphological  transition  in  the 
subaerial  and  subaqueous  nearshore  zone  at  Duck,  North 
Carolina  is  studied  in  this  report.  Specifically,  we 
examine  the  role  of  various  synoptic  weather  systems  as 
agents  of  change  in  the  sedimentary  prism;  the  importance  of 
wave  height,  period,  direction,  steepness  and  breaker  type 
relative  to  profile  changes;  and  the  role  played  by  existing 
profile  form  on  this  seasonal  change  in  nearshore  mor¬ 
phology. 


LITERATURE  REVIEW 

The  literature  on  causes  of  changes  in  beach  profiles 
can  be  divided  into  contributions  that  focus  on  the  role  of 


pre-existing  morphology,  storms,  incident  waves,  and  edge 
waves.  The  concept  that  morphological  changes  are  a 
function  of  the  pre-existing  morphology  implies  an  ordered 
sequence  of  change.  Sonu  and  Van  Beek  (1971)  offered  such  a 
model  of  sequential  beach  changes  based  on  data  collected 


along  the  Outer  Banks  o£  North  Carolina  (Dolan,  1966),  and 
they  proposed  a  linkage  between  rhythmic  topography  and 
shore-normal  profile  change.  Their  along-the-coast  study 
reach  was  limited  to  tens  of  meters,  however.  The  model 
they  developed  allowed  the  prediction  of  successive  beach 
profiles  based  on  pre-existing  sediment  storage,  beach  width 
and  surface  configuration. 

More  recently  Wright  et  al.  (1979)  proposed  a  model  of 
sequential  beach  changes  relating  to  surf  zone  morphody- 
namics  based  on  studies  of  Australian  beaches.  Two  major 
types  of  beaches  were  identified:  steep,  reflective  beaches 
and  flat,  dissipative  beaches.  Steep,  reflective  beaches 
have  well-developed  berms  and  beach  cusps,  and  surging 
breakers  with  high  runup  and  minimum  setup.  Rip  currents  and 
the  associated  three  dimensional  inshore  topography  are 
absent  in  steep,  reflective  beaches.  Flat,  dissipative 
beaches  have  wide  surf  zones,  much  more  complex  nearshore 
topography  (than  reflective  systems),  and  rip  cells.  Short 
(1979),  also  working  on  Australian  beaches,  observed  sequen¬ 
tial  bar  stage  development  which  he  described  in  a  three- 
dimensional  model  of  ten  stages.  These  changes  in  the  ob¬ 
served  beach  morphology  were  related  to  the  cumulative  wave 
power . 

Storm  waves  have  been  recognized  as  modifying  agents 
on  beaches  for  many  years  (Zeigler,  et  al.,  1959;  Fox  and 
Davis,  1976;  Dolan  and  Hayden,  19» l) .  Zeigler  et  al.  (1959) 


found  that  storms  along  Outer  Cape  Codr  Massachusetts, 
caused  rapid  changes  in  beach  volume  and  topography.  Their 
study  concluded  that  after  storms,  profiles  were  either 
planar  or  concave  upward,  while  during  quiescent  periods 
they  were  convex  upward.  In  a  more  recent  study,  Fox  and 
Davis  (1976)  investigated  the  responses  of  coastal  systems 
to  atmospheric  (storm)  processes.  Working  with  data  from 
the  eastern  shore  of  Lake  Michigan;  Mustang  Island,  Texas; 
and  the  central  Oregon  coast,  they  developed  a  coastal 
process-response  model.  Time  series  analyses  of  weather, 
waves  and  current  data  were  combined  with  frequent  beach 
profile  data  to  construct  the  model.  The  study  concluded 
that  morphologic  responses  differed  in  the  three  areas 
because  their  positions  relative  to  storm  tracks,  coastal 
orientation,  tidal  range,  and  nearshore  topography  were 
d'  ferent.  Dolan  and  Hayden  (1981)  examined  the  along- 
the-coast  periodicities  in  storm  overwash  following  the  Ash 
Wednesday  Storm  of  March  1962.  This  study  concluded  that 
the  location  and  magnitude  of  storm  deposits  and  storm 
hazards  were  systematically  distributed  along  the  Atlantic 
Coast. 

Goldsmith  et  al.  (1982)  delineated  sequential  bar 
development  patterns  at  Hahoterim  Beach,  Israel.  These 
changes  were  found  to  be  induced  by  changes  in  the  incident 
wave  energy  of  the  area,  and  the  threshold  values  for  these 
morphologic  transitions  were  quantified.  They  concluded 


that  double-crescentic  patterns  were  the  optimal  evolution¬ 
ary  stage  of  crescentic  systems,  and  that  outer  bars  and 
mega-rip  channels  composed  the  most  stable  inshore  elements. 
Aubrey  et  al.  (1980)  developed  a  linear  statistical  estima¬ 
tion  model  (based  on  empirical  eigenfunctions)  to  predict 
changes  in  beach  profiles,  using  variations  in  incident  wave 
energy.  They  concluded  that  weekly  mean  wave  energy  was  the 
best  predictor  of  beach  changes. 

Rhythmic  coastal  topography  has  been  explained  in 
several  ways:  intersecting  wave  trains,  variations  in 
nearshore  topography,  and  edge  waves.  The  concept  of  edge 
waves  as  morphologic  agents  has  been  investigated  by  many 
researchers  (Dalrymple,  1975;  Bowen  and  Inman,  1971;  Guza 
and  Inman,  1975;  Dolan  et  al.,  1979).  Edge  waves  are  normal 
trapped  modes  of  longshore  periodic  wave  motions.  These 
motions  occur  along  the  edge  of  water  bodies  and  may  be  one 
of  two  distinct  types:  standing  or  progressive.  These 
phenomena  can  be  found  on  straight  or  curving  shorelines,  as 
well  as  on  concave  or  convex  offshore  slopes.  Their 
inability  to  radiate  energy  into  deeper  water  causes  them  to 
dissipate  their  energy  only  through  friction  and  interaction 
with  other  currents  and  wave  forms  (Guza  and  Inman,  1975) . 

Guza  and  Inman  (1975)  analyzed  beach  cusps  in  reflec¬ 
tive  systems  from  both  analytical  and  observational  ap¬ 
proaches.  They  found  that  the  edge  wave  resonances, 
theoretically  predicted  and  observed,  were  not  visible  on 
the  beach  face  during  plunging  incident  wave  conditions. 


In  a  more  recent  study,  Dolan  et  al.  (1979)  used  the 
edge  wave  phenomenon  to  explain  the  longshore  periodicities 
along  the  Cape  Hatteras  Arc.  They  noted  that  standing  edge 
waves,  unlike  intersecting  wave  trains  and  inshore  topogra¬ 
phy,  would  result  in  both  temporal  and  spatial  stability  of 
the  shore  face  features.  They  concluded  that  edge  waves 
were  important  in  determining  the  intensity  of  regional- 
scale  morphodynamic  responses  along  the  mid-Atlantic  coast 
barrier  islands. 


STUD*  SITE 

The  Coastal  Engineering  Research  Center  (CERC)  has 
been  conducting  coastal  research  for  over  twenty  years. 
Field  observations  supportive  of  earlier  experimental  and 
theoretical  investigations  lacked  accurate,  repetitive  wave, 
beach  and  water  level  data,  however,  especially  during 
storms.  To  meet  the  need  for  such  data  the  Field  Research 
Facility  (FRF)  was  constructed  near  Duck,  North  Carolina 
(Fig.  1).  Located  at  36°10 '  north  latitude  and  75°45'  west 
longitude,  the  FRF 1 s  561  m  research  pier  and  the  418  m2 
laboratory  and  office  building  were  completed  in  1976  and 
1980,  respectively  (Birkemeier,  et  al.,  1981).  The  addition 
of  the  Coastal  Research  Amphibious  Buggy  (CRAB)  to  the  FRF's 
data  collection  facilities  in  1981  has  improved  the  quality, 
density  and  resolution  of  CERC’s  profile  data  base. 

The  coast  at  the  FRF  is  part  of  an  unbroken  100  km 
stretch  of  coastline  extending  south  from  Rudee  Inlet, 


Point  Harbor 


Virginia,  to  Oregon  Inlet,  North  Carolina.  This  long, 
narrow  barrier  island  is  part  of  the  Outer  Banks,  North 
Carolina.  The  area  is  subject  to  winds  and  waves  from 
intense  high  pressure  cells  as  well  as  tropical  and  extra- 
tropical  low  pressure  cells.  The  climate  is  temperate 
marine  with  mild  winters  and  warm  summers  (Birkemeier  et 
al.,  1981). 

Forewarning  of  storms  is  difficult  in  this  area 
because  cyclogenesis  often  occurs  offshore  from  Cape  Hat- 
teras.  North  Carolina.  Bosserman  and  Dolan  (1968)  deter¬ 
mined  that  the  most  damaging  storms  follow  northeast  tracks 
off  the  mid-Atlantic  coast.  The  highest  winds  and  waves 
tend  to  occur  when  low  pressure  storm  movement  is  blocked  by 
a  high  pressure  system  to  the  north.  Hayden  (1981) 
calculated  the  mean  extratropical  cyclone  frequency  for 
Eastern  North  America  in  2.5°  latitude  by  5.0°  longitude 
grid  cells.  His  results  showed  a  storm  frequency  maximum 
along  the  east  coast  baroclinic  zone.  Tropical  cyclone 
frequencies  along  the  North  Carolina  coast  are  variable. 
The  area  between  Cape  Lookout  and  Cape  Hatteras  has  the 
highest  frequency  of  hurricane  occurrence,  while  the  area 
around  the  FRF  has  the  lowest. 

STORM-WAVE  EVENTS 

During  the  nine-month  survey  period,  six  synoptic 
scale  weather  events  occurred  that  caused  winds  of  suffi¬ 
cient  magnitude  to  generate  wave  heights  in  excess  of  2  m  at 


the  Field  Research  Facility.  These  events  were  divided  into 
three  major  categories:  tropical  storms/hurricanesf  extra- 
tropical  storms,  and  high  pressure  systems.  There  were  two 
examples  of  each  storm  type  during  the  study  period.  The 
significant  environmental  variables  for  each  of  the  events 
are  summarized  in  Table  2. 

Hurricane  Dennis,  the  first  1981  tropical  storm  to 
effect  the  study  area,  originated  from  an  African  distur¬ 
bance  and  was  named  on  August  8,  1981,  while  centered  800  km 
southwest  of  the  Cape  Verde  Islands.  The  system  tracked 
westward  across  the  Atlantic  and  Caribbean,  nearly  stalling 
on  August  14  south  of  Cuba.  On  August  15,  Dennis  reinten¬ 
sified  and  turned  northward  across  Cuba  and  toward  the 
southeastern  United  States.  Dennis  travelled  parallel  to 
the  east  coast  before  turning  northeastward  away  from  the 
middle  Atlantic  states  on  August  20.  The  highest  sustained 
surface  wind  reported  by  a  land  station  was  20  m/s  at  Cedar 
Island,  North  Carolina.  Dennis  reached  hurricane  status  for 
approximately  12  hours  on  August  20,  after  which  the  storm 
center  moved  over  colder  waters  and  lost  its  tropical 
characteristics  by  August  22.  Maximum  sustained  winds  were 
36  m/s  and  the  lowest  central  pressure  was  995  mb  (Lawrence 
and  Pelissier,  1982) 

Hurricane  Emily  formed  on  a  frontal  wave  in  a  manner 
similar  to  Cindy  (August  2-5,  1981).  The  storm  tracked 
northeastward  toward  Bermuda.  The  storm  system's  forward 


progress  was  temporarily  impeded  by  a  high  pressure  system 
passing  to  the  north,  causing  the  tropical  storm  to  make  a 
small,  counterclockwise  loop  on  September  3,  1981.  It  was 
during  this  loop  phase  that  Emily  tracked  westward  and  posed 
its  greatest  threat  to  the  middle  Atlantic  region.  On 
September  4,  Emily  began  to  move  northeastward  and  back  out 
to  sea.  Maximum  surface  winds  of  41  m/s  (and  a  minimum 
surface  pressure  of  966  mb)  were  reached  on  September  6, 
(Lawrence  and  Pelissier,  1982). 

On  October  11,  1981,  the  first  of  two  consecutive 
strong  high  pressure  systems  formed  in  the  Maritime  Prov¬ 
inces  of  Canada.  This  high  intensified  and  moved  into  the 
New  England  States  on  the  13th.  A  strong  pressure  gradient 
and  onshore  winds  prevailed  at  the  PRF  throughout  the 
duration  of  the  event.  By  the  15th,  the  ridge  had  weakened 
and  the  gradient  deteriorated  as  the  high  moved  off  the 
coast  of  Maine.  The  maximum  pressure  reading  (1036  mb)  was 
observed  on  October  13,  at  1200  GMT. 

A  second,  stronger  high  pressure  cell  formed  in  the 
same  region  as  the  first  on  October  29,  1981.  This  system 
also  produced  strong  pressure  gradients  normal  to  the  east 
coast,  with  persistent  onshore  winds.  The  maximum  pressure 
(1042  mb)  was  observed  on  October  31,  while  the  high  was 
centered  over  northern  Maine.  The  system  dissipated  as  it 
tracked  eastward  off  the  New  England  coast. 

The  first  extratropical  storm  event  of  the  season 
formed  over  the  southwestern  United  States.  The  storm 


travelled  southeastward  Into  the  Gulf  of  Mexico  and  tra¬ 
versed  Florida  on  November  11,  1981.  By  November  13,  the 
storm  was  centered  off  Cape  Hatteras  and  ships  in  this  area 
were  reporting  40-50  knot  winds  with  waves  in  excess  of  30 
feet.  At  7:00  A.M.  EST  on  the  15th,  the  storm  was  centered 
at  approximately  38°N  latitude  and  72°W  longitude  and  had  a 
central  pressure  of  984  mb.  The  combination  of  the  low 
pressure  cell  and  a  syzygy-perigean  alignment  of  the  sun, 
moon  and  earth  resulted  in  high  waves  and  water  levels  at 
the  study  site  (Miller  and  Leffler,  1981).  The  storm 
continued  to  track  northeastward,  crossing  Newfoundland  by 
the  18th  (Mariner's  Weather  Log,  1982). 

On  November  24,  1981,  a  secondary  low  was  triggered 
off  the  South  Carolina  coast  as  a  low  pressure  cell  moved 
into  the  Appalachian  Mountains  of  the  Carolinas.  By  7:00 
A.M.  on  the  25th,  the  storm  had  depressed  to  980  mb  and  was 
located  off  Cape  Hatteras.  An  offshore  buoy  recorded  26 
foot  wave  heights,  and  winds  at  Cape  Hatteras  exceeded  50 
knots.  By  November  26,  the  storm  was  972  mb,  but  it  was 
located  well  out  to  sea  (40°N,  57°W)  (Mariner's  Weather  Log, 
1982) . 

DATA  SET 

Profile  data  were  collected  at  the  Field  Research 
Facility  (FRF),  Duck,  North  Carolina,  by  the  O.S.  Army  Corps 
of  Engineers,  using  the  Coastal  Research  Amphibious  Buggy 


(CRAB).  The  CRAB  is  a  10.7  m  high  amphibious  tripod  vehicle 
capable  of  operating  in  water  depths  up  to  9  m  and  wave 
heights  up  to  1.8  m.  Position  and  profile  elevation  are 
determined  using  the  CRAB-Zeiss  surveying  system;  a  Zeiss 
Elta-2  first  order,  self-recording  electronic  distance  meter 
in  combination  with  the  CRAB  (Birkemeier  et  al.,  1981  and 
Miller  and  Leffler,  1982). 

Ten  bathymetric  surveys  containing  11  individual  pro¬ 
file  lines  were  evaluated  in  this  study  (Fig.  2).  The  study 
period  spanned  nine  months  between  July  2,  1981,  and  March 
18,  1982  (Table  3).  Environmental  variables  for  the 
corresponding  storm  periods  were  also  collected  at  the  FRF, 
and  reported  in  the  monthly  Basic  Environmental  Data 
Summary,  published  by  CERC.  Daily  values  for  barometric 
pressure,  wind  speed,  wind  direction,  wave  height,  wave 
period,  and  longshore  current  direction  were  used  in  this 
study. 

The  profile  data  were  organized  for  analysis  in  two 
ways:  shore  parallel  transects  and  shore  normal  transects 
(traditional  method).  The  profile  lines  normal  to  the  shore 
each  contain  44  elevation  points  (variables)  beginning  60  m 
seaward  of  the  FRF  baseline  and  extending  to  760  m  from  the 
baseline.  This  entire  data  matrix  is  44  variables  by  110 
cases.  The  shore  parallel  data  are  structured  as  an  11 
variable  by  440  case  data  matrix.  The  shore-parallel  data 
matrix  is  further  subdivided  into  two  groups;  one  for 
profiles  north  of  the  pier  and  the  other  for  profiles  south 
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July  17,  1981 
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August  4,  1981 

Survey 

4s 

August  24,  1981 

Survey 

5* 

September  19.  1981 

Survey 

6: 

November  3.  1981 

Survey 

7t 

November  16 ,  I98I 

Survey 

8: 

January  5.  1982 

Survey 

9« 

February  9,  1982 

Survey 

10  s 

March  18,  1982 

TABLE  3 :  Bathymetric 
the  FRF. 
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of  the  pier.  These  profile  data  sub-sets  contain  220  cases 
each  and  six  and  five  variables,  respectively. 


THE  ANALYSIS 

Principal  components  analysis  was  used  to  define  the 
modes  of  intercorrelation  that  exist  among  the  variables 
(profile  elevation  points)  .  The  analysis  was  designed  to 
delineate  both  the  shore-parallel  and  the  shore-normal 
variations  in  nearshore  morphology. 

Principal  components  analysis  has  been  successfully 
used  to  partition  the  variance  in  multivariate  geophysical 
data  sets  (Kutzbach,  1967;  Resio  et  al.,  1974;  Hayden  et 
al.,  1975;  Aubrey,  1979;  Aubrey  et  al.,  1980;  Hayden,  1981; 
Dolan  et  al.,  1982).  Using  a  minimization  of  least  square 
errors,  this  method  of  analysis  is  useful  for  determining 
patterns  in  large  data  sets  (Lorenz,  1956;  Gilman,  1957; 
Kutzbach,  1967).  Principal  components  analysis  transforms  a 
series  of  intercor related  variables  into  a  set  of  new 
statistically  independent  variables.  These  new  variables 
are  linear  combinations  of  the  original  variables  but  are 
mutually  orthogonal.  Principal  components  analysis  was 
chosen  in  this  study  because  profile  elevation  variations 
are  quantified  and  organized  in  the  first  few  principal 
components,  and  because  the  principal  component  eigenvectors 
and  their  weightings  in  time  can  be  related  directly  to  the 
physical  environment. 
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RESULTS 


Means  and  Standard  Deviations 

Fig.  3  shows  the  mean  profile  elevations  along  with 
standard  deviations  of  these  elevations  at  each  of  the  44 
intervals  along  the  profile  line.  A  single  bar  is  located 
approximately  340  m  from  the  baseline  (station  number  27). 
The  remainder  of  the  profile  is  featureless.  Maximum 
standard  deviations  occur  in  the  center  of  the  profile 
(points  17-21) ,  and  decrease  both  landward  and  seaward.  As 
noted  in  the  Horn  Island,  Mississippi  data  collected  by 
Dolan  et  al.  (1982),  there  is  no  variance  maxima  in  the 
vicinity  of  the  bar  (point  27).  The  means  and  standard 
deviations  of  the  profile  elevation  points  are  given  in 
Table  4. 

Fig.  4  shows  the  means  of  the  11  profiles  at  each 
survey  date.  At  the  time  of  the  July  2,  1981,  survey,  a 
single  inshore  bar  was  centered  220  m  from  the  baseline. 
The  July  17,  1981,  and  August  4,  1981,  mean  profiles  showed 
a  slight  filling  of  the  inner  trough.  Hurricane  Dennis' 
passage  caused  a  seaward  displacement  of  the  bar  out  to  280 
m  from  the  baseline  in  the  August  24,  1981.  survey.  The 
trough  also  deepened  slightly  during  this  time. 

By  the  September  19,  1981,  survey,  a  double-bar 
profile  had  evolved  with  a  broad  inner  bar  centered  100  m 
from  the  baseline  and  a  small  outer  bar  located  approxi¬ 
mately  300  m  from  the  baseline.  There  was  a  substantial 
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FIGURE  3  :  The  mean  profile  plus-or-minus  one  standard  deviation  at  the  FRF 


th  c\j  \t\\q  o-co  on  o  t-i  cj  u%vo  c^oo  on  o  »h  cm  rv^-  vr>vo  o-oo  on  o  ih  cm  rvt  nano  o-co  on  o  »-t  cnj  c»v}- 

THTHrHrMTH'HT-l'H.H'rMCMOJCMOJCMCNjCNJCNICNJCSJ  CA  CA  CA  CA  CA  (A  CA  CA  CA-^ ^ 


OFFSHORE 

STATION 

NUMBER 


DISTANCE 

FROM 

BASELINE* 


DISTANCE 

OFFSHORE* 


MEAN* 


STANDARD 

DEVIATION 


*  These  values  are  expressed  in  meters  (m) . 

TABLE  4  :  The  mean  elevation  (measured  from  MSL)  and  standard 
deviation  for  each  of  the  44  stations  along  shore- 
normal  profiles.  Station  locations  are  measured 
from  both  the  CERC  baseline  and  the  shoreline. 
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FIGURE  4  :  The  mean  FRF  profiles  for  each  survey  date.  The 
profiles  are  offset  by  1  (one)  meter  intervals. 


amount  of  erosion  in  the  berm  area  following  this  survey. 
The  November  3,  1981,  mean  profile  had  a  well  developed 
storm  bar  300  m  from  the  baseline,  with  a  broad  inner 
trough.  Berm  destruction  continued  through  this  survey 
period.  Maximum  development  of  the  storm  bar  occurred  by 
the  November  16.  1981,  survey.  The  outer  bar  had  moved  60  m 
farther  offshore  and  a  small  inner  bar  had  developed  140  m 
from  the  baseline. 

The  survey  of  January  5,  1982,  showed  an  increase  in 
size  of  the  inner  and  outer  bar  systems,  and  both  bars  were 
slightly  farther  offshore  relative  to  the  previous  survey. 
The  February  9,  1982,  mean  profile  showed  no  change  seaward 
of  the  outer  bar,  while  sediment  deposition  occurred  land¬ 
ward  of  this  point.  The  final  survey,  March  18,  1982, 
showed  a  slight  broadening  of  the  outer  bar  with  little  or 
no  change  in  the  nearshore  areas. 

The  Eigenvectors 

The  percent  of  variance  explained  in  each  of  the  first 
four  eigenvectors  of  the  correlation  matrix  is  summarized  in 
Table  5  for  shore  normal  profile  data.  The  Monte  Carlo  test 
of  Overland  and  Preisendorfer  <1982)  for  statistical  signif¬ 
icance  of  the  eigenvectors  was  used  here.  The  first  four 
eigenvectors  explained  84.4%  of  the  total  variance.  Each  of 
these  vectors  explained  a  greater  percentage  of  the  variance 
than  are  explained  by  95%  of  the  same  order  vectors  computed 
from  random  data  matrices  of  the  same  dimensions.  Geophys- 


ical  significance  was  determined  by  inspection  and  study  of 
each  individual  eigenvector  and  its  variation  in  space  and 
time. 

The  first  four  eigenvectors,  those  that  are  statis¬ 
tically  significant  at  the  .05  level,  were  examined  (Fig. 
5a,  6a.  7a,  and  8a).  To  aid  in  the  interpretation  of  the 
eigenvectors,  profile  elevations  were  reconstructed  by  ad¬ 
ding  the  variation  explained  by  each  eigenvector  according 
to 

j  *  X^  ±C  (SDj)  (E^ j ) 

where  is  the  reconstructed  value  for  the  ith  point  along 
the  profile  of  the  jth  eigenvector.  X i  is  the  mean  profile 
elevation  at  the  ith  point  along  the  profile.  SD^  is  the 
standard  deviation  of  profile  elevation  at  the  ith  point 
along  the  profile,  and  E^j  is  the  weighting  on  the  jth 
eigenvector  for  the  ith  point  along  the  profile.  ±C  is  a 
constant  which  is  representative  of  the  case  weightings 
(time  series)  on  a  particular  vector.  Plots  of  R^j  at  each 
profile  point  for  representative  positive  and  negative  C 
values  (Table  6)  are  shown  in  Fig.  5b,  6b,  7b,  and  8b. 

The  average  case  loading  (Afc)  for  the  jth  eigenvector 
at  profile  location  over  all  survey  dates  is  given  by 

1  10 

At“  io  z  E^k 

t-1 

at  each  k  and  for  j«l  to  4.  These  case  loadings  are  shown  in 
Fig.  5c,  6c,  7c,  8c. 
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PROFILE  LINE  NUMBER  SURVEY  NUMBER 

FIGURES  8a  -  8d:  a)  Plotted  values  for  the  fourth  vector  (E4) 5  b)  Reconstructed  pro 

file  for  E4  plotted  plus  the  mean  ( - )  and  minus  the  mean  (...); 

c)  Mean  vector  weighting  on  E4  over  all  survey  dates  for  each  in¬ 
dividual  profile;  d)  Mean  vector  weighting  on  E4  for  all  profiles 
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The  average  case  loading  Ak  for  the  j  eigenvector  at 
the  tfck  CRAB  survey  over  all  profile  locations  (k)  is  given 
by 

1  11 

A* '  n  1 
k»l 

at  each  survey  date  t,  for  j  *  1  to  4.  These  case  loadings 
are  shown  in  Fig.  5d,  6d,  7d,  and  8d. 

The  first  eigenvector  (El)  explained  35.4%  of  the 
total  variance  and  represented  the  dominant  mode  of  profile 
variations  from  the  mean.  A  positive  case  weighting  on  El 
indicates  a  deficiency  of  sediment  on  the  beach  face  out  to 
300  m  from  the  baseline  (Fig.  5a).  Seaward  of  this  point, 
there  is  a  surplus  of  sediment.  A  negative  case  weighting 
on  El  indicates  more  sediment  in  the  inshore  region  and  less 
in  the  offshore  area.  El  is  a  profile  rotation  function 
with  positive  weightings  when  the  profile  is  flatter  than 
the  mean,  and  negative  when  steeper.  While  it  was  tempting 
to  attribute  the  along- the-shor e  variation  in  mean  El 
weightings  (Fig.  5c)  to  the  presence  of  the  pier  (Miller  et 
al.,  1983),  the  along-the-coast  distance  between  the  maximum 
weighting  at  profile  line  number  2  to  the  minimum  at  number 
7  was  343  m?  and  Vincent  (1973)  showed  that  crescentic  sand 
waves  with  wave  lengths  of  140  to  650  m  are  common  along  the 
North  Carolina  coast.  The  meanders  in  the  crescentic  bar 
forms  were  slightly  skewed  to  the  north  as  were  the  mean 


loadings  on  El  calculated  in  this  study.  This  study  did  not 
affix  cause,  only  association;  and  it  is  clear  that 
along-the-coast  profile  slope  variations  at  this  length 
scale  do  occur. 

The  time  history  of  El  case  loadings  for  the  average 
profile  location  is  shown  in  Fig.  5d.  Positive  values 
occurred  from  the  first  through  the  fifth  survey.  After  the 
fifth  survey  (September  19,  1981),  negative  weightings 
prevailed.  The  transition  from  a  relatively  steep  to  a 
flatter  profile  slope  was  abrupt. 

The  second  eigenvector  (E2)  (Fig.  6a)  explained  25.7% 
of  the  total  variance.  A  positive  weighting  on  E2  indicates 
a  single  offshore  bar  at  250  m  (Fig.  6b).  A  negative 
weighting  indicates  a  seaward  displacement  of  sediment,  the 
development  of  an  inner  bar  140  m  from  the  baseline,  and  an 
outer  bar  at  360  m. 

The  maximum  average  loading  on  E2  over  all  survey  data 
was  located  on  profile  line  number  1,  north  of  the  pier 
(Fig.  6c).  Average  loadings  declined  from  this  profile 
line,  then  increased  again  on  the  south  side  of  the  pier. 
There  was  a  tendency  for  the  single-bar  configuration  to 
occur  at  the  northern  and  southern  ends  of  the  study  site 
and  for  the  double-bar  configuration  to  occur  in  the  center 
of  the  reach.  The  spatial  scale  of  this  variation 
suggested  an  along-the-coast  variability  greater  than  sev¬ 
eral  km. 


The  temporal  history  of  the  vector  loadings  (averaged 
over  all  profile  locations)  (Fig.  6d)  indicated  a  change 
from  the  single-bar  system  prior  to  the  6th  survey  (November 
3rd)  to  the  double-bar  configuration  thereafter.  As  was  the 
case  with  El*  the  transition  in  temporal  loadings  on  E2  was 
abrupt  and  occurred  between  September  and  November. 

The  third  eigenvector  (E3)  (Fig.  7a)  accounted  for 
13.1%  of  the  total  variance.  A  positive  weighting  on  E3 
indicates  a  bar  290  m  from  the  baseline*  a  broad,  deep  inner 
trough  and  less  sediment  on  the  beach  face  (Fig.  7b).  A 
negative  weighting  indicates  buildup  of  the  beach  face*  a 
broad  inner  bar,  and  a  surplus  of  sand  seaward  of  360  m. 


The  mean  loading  on  each  profile  for  E3  was  positive, 
except  on  three  profiles  to  the  south  of  the  pier  (Fig.  7c). 
The  time  history  of  E3  was  marked  by  a  period  of  positive 
loadings  which  occurred  simultaneously  with  the  changes  in 
El  and  E2  (Fig.  7d) . 

Eigenvector  four  (E4)  (Fig.  8a)  explained  10.2%  of  the 
total  variance.  A  positive  weighting  on  E4  indicates  an 
inner  bar  200  m  from  the  baseline  and  a  deficiency  of  sand 
along  the  remainder  of  the  profile,  except  for  the  extreme 
offshore  region  seaward  of  500  m  (Fig.  8b).  A  negative 
weighting  indicates  a  sand  surplus  at  the  beach  face,  a 
well-developed  outer  bar  located  approximately  300  m  from 
the  baseline,  and  a  trough  in  between. 
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The  spatial  variations  in  E4  were  small  (Fig.  8c). 
The  majority  of  the  profiles  to  the  north  of  the  pier  were 
positively  weighted,  whereas  the  southern  profiles  were 
negatively  weighted.  The  time  series  of  the  E4  loadings 
(Fig.  8d)  were  very  significant.  The  weightings  were 
positive  on  E4  through  the  September  19,  1981,  survey,  after 
which  the  weightings  became  negative.  This  coincided  with 
the  major  shifts  in  the  other  three  vectors. 

Storms  and  Profile  Changes 

During  the  study  period  three  different  storm  types 
occurred:  hurricanes/tropical  storms,  extratropical  cyclones 
and  anticyclones.  Table  2  summarizes  the  important  parame¬ 
ters  of  the  six  individual  storm  events. 

Hurricane  Dennis  and  Tropical  Storm  Emily  caused  only 
minor  changes  in  the  FRF  bathymetry.  Dennis  caused  varia¬ 
tion  only  on  profiles  4  and  6  (as  is  evidenced  in  the  E1-E2 
plots.  Fig.  9-10).  The  effects  on  the  two  profiles  were 
opposite:  positive  on  4  and  negative  on  6.  Tropical  Storm 
Emily  caused  sufficient  changes  on  several  profiles  to 
significantly  alter  the  mean  profile  for  the  September  19, 
1981,  survey  (Fig. 4).  The  beginning  of  a  small  offshore  bar 
appeared  in  this  survey. 

Two  high  pressure  systems  occurred  in  October,  1981. 
generating  high  winds  and  waves  at  the  pier.  During  this 
month,  large  changes  were  made  in  the  local  bathymetry.  The 


FIGURES  10a  -  10b:  Profile  lines  7  through  11  plotted  in  E1-E2  eigenvector  space 


largest  changes  in  E3  and  E4  were  also  observed  during  this 
period,  along  with  the  second  largest  variations  in  El  and 


E2.  The  morphological  responses  to  these  changes  were 
offshore  sediment  movement  and  the  development  of  the  storm 
bar  approximately  300  m  from  the  FRF  baseline. 

During  November,  the  two  extratropical  cyclones  that 
skirted  the  mid-Atlantic  coast  made  additional  changes  in 
the  profiles.  The  first  storm  caused  large  changes  in  El 
and  E2  and  smaller  changes  in  E3  and  E4.  This  storm  moved 
the  storm  bar  60  m  farther  offshore  and  further  developed 
the  inner  bar.  The  second  storm  continued  the  development 
of  storm  profile  features. 

Profile  Form  and  Profile  Changes 

Two  survey  periods  (September  19,  1981,  to  November  3, 
1981;  and  November  3  to  November  16)  were  analyzed  to  answer 
questions  about  the  importance  of  original  profile  form  on 
subsequent  profile  change.  The  bulk  of  the  profile  changes 
during  the  study  period  occurred  during  these  surveys. 
These  changes  were  analyzed  by  plotting  the  loading  on  a 
vector  at  one  survey  date  against  the  change  in  the  vector 
loading  at  the  next  survey.  The  graphs  (Fig.  11  and  12) 
indicated  onshore-offshore  movement  of  sediment  (El  vs  A  El) 
and  the  development  of  bars  on  the  profile  (E2  vs  AE2) . 
Linear  regression  was  used  to  assign  a  degree  of  correlation 
and  significance  to  the  relationship. 


1(9/19  to  11/3)  |  |  J_aE1(11/3  to  11/16) 
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Fig.  11a  shows  the  changes  in  El  during  October  as  a 
function  of  El  at  the  beginning  of  the  period.  This  was  the 
period  in  which  two  high  pressure  systems  were  responsible 
for  two  large  wave  events.  The  mean  profiles  showed 
positive  changes  in  El  for  increasing  values  of  El  at  the 
beginning  of  the  period.  A  positive  change  in  the  weight¬ 
ings  on  El  indicates  offshore  sediment  movement,  while  a 
negative  change  indicates  onshore  sediment  movement.  The 
north-of-pier  profiles  had  a  (El  vs  AE1)  correlation 
coefficient  of  .593  <a95  *  .11) .  Thus,  the  flatter  the 
profile,  the  greater  is  the  subsequent  offshore  sediment 
movement,  further  flattening  the  profile.  The  changes  in  El 
for  the  period  spanning  the  occurrence  of  the  two  cyclones 
of  November  are  shown  in  Fig.  lib  relative  to  the  weighting 
on  El  at  the  beginning  of  the  period.  The  correlations 
between  El  and  AEl  for  the  areas  to  the  north  and  south  of 
the  pier  were  -.965  <<*95*. 0009)  and  +.832,  (a95  •  .04), 
respectively.  South  of  the  pier,  offshore  displacement  of 
sediment  was  greater  in  areas  with  flatter  profiles  prior  to 
the  survey  period.  The  opposite  relationship  prevailed 
north  of  the  pier. 

Fig.  12  shows  the  changes  in  E2  relative  to  E2  at  the 
beginning  of  the  period  of  high  pressure  systems  in  October. 
Negative  weightings  on  E2  indicate  the  tendency  for  a 
double-bar  profile.  Profiles  with  a  more  highly  developed 
original  single-bar  configuration  showed  a  greater  tendency 
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to  change  toward  a  double-bar  configuration.  The  north- 
of-pier  data  points  had  a  correlation  coefficient  of  -.760, 
and  the  combined  data  points  showed  a  correlation  of  -.650 
(095**02).  The  changes  in  E2  during  the  passage  of  the 
extratropical  cyclone  in  November  showed  negative  displace¬ 
ment,  indicating  storm  profile  development  (Fig.  12b).  The 
changes  were  greater  for  this  latter  period  than  in  the 
earlier  period.  The  south-of-pier  data  points  showed  a 
correlation  coefficient  of  -.641  (a^5«.12).  From  this  we 
inferred  that  for  both  the  October  high  pressure  period  and 
the  November  cyclone  period,  the  change  toward  two  bars  was 
most  extreme  where  the  single-bar  form  was  most  highly 
developed. 

The  relationship  between  E2  for  a  survey  and  the 
subsequent  change  in  E2  was  analyzed  for  each  survey 
interval  (Fig.  13  and  Table  7).  A  statistically  significant 
relationship  was  found  for  all  those  intervals  in  which 
storms  occurred.  In  each  case  the  slope  of  the  relationship 
was  negative.  Apparently  the  slope  is  a  reflection  that 
incident  wave  energy  from  storms  was  applied  in  these 
periods,  work  was  done  on  the  sedimentary  prism,  and  large 
E2s  resulted.  The  x-intercept  indicates  the  equilibrium  E2 
weighting  for  the  period.  For  the  period  that  included 
Hurricane  Dennis  the  equilibrium  value  of  E2  was  3.3. 
Profiles  prior  to  Dennis  with  E2s  greater  than  3.3  had  a 
decline  in  E2  during  the  period  while  those  profiles  with  an 


FIGURE  13:  Regression  lines  for  correlations  between  the  weightings  on  E2  and  the  change 

in  E2  between  survey  periods. 
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E2  less  than  3.3  had  an  increase  in  E2  over  the  period.  For 
the  survey  spanning  tropical  storm  Emily  the  equilibrium  E2 
weighting  was  3.7.  For  the  period  of  the  two  high  pressure 
systems  in  October  the  equilibrium  condition  was  0.8, 
similar  to  the  grand  mean  profile  but  with  a  greater 
amplitude  in  the  single-bar-and-trough  configuration.  In 
the  November  storm  period  the  equilibrium  point  for  E2  was 
-6.4.  The  two-bar  profile  was  the  equilibrium  condition. 
The  current  study  was  not  designed  to  answer  the  question  of 
the  cause  of  a  particular  equilibrium  "setting”  for  E2,  but 
we  will  offer  some  speculation  in  the  discussion  section. 

Profile  Changes  In  E1-E2  and  E3-E4  Space 

At  each  survey  date,  each  profile  may  be  represented 
as  a  point  in  eigenvector  space.  The  results  of  the 
analysis  in  E1-E2  space  are  shown  in  Table  8  and  Fig.  9,  10 
and  14.  For  E3-E4  space  the  results  are  shown  in  Fig.  15 
and  16.  A  plot  of  the  temporal  weightings  on  El  and  E2  (E3 
and  E4)  in  a  cartesian  coordinate  system  shows  the  nature  of 
change  in  profile  shape  between  survey  dates.  The  length  of 
the  individual  line  segments  connecting  these  points  is  a 
measure  of  the  magnitude  of  change  in  profile  morphology 
between  surveys.  By  these  plots  the  periods  of  maximum 
profile  change  at  the  FRF  can  be  easily  identified.  To 
better  understand  these  sediment  movements  in  eigenvector 
space,  the  characteristic  profile  for  the  four  cartesian 


«  July  2,  1981 


TABLE  8 


i  July  17,  198I 

*  August  4,  1981 
»  August  24,  198I 

*  September  19,  1981 

:  November  3.  1981 

1  November  16,  1981 

1  January  5,  1982 

1  February  9,  1982 

1  March  18,  1982 


s  Summary  of  symbols  used  in  FIGURES  9,  10,  and 
15a  through  15k. 
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FIGURE  l4:  Characteristic  profiles  in  each  of  the  cartesian  quadrants  for  the  plots 
in  E1-E2  eigenvector  space. 


FIGURES  15a  -  15d:  Profile  lines  1  through  4  plotted  in  E3-E4 

eigenvector  space. 
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quadrants  is  reconstructed  as  follows: 

Rij  -  Xi  +  (iCHSDiXEij)  +  (±C)  (SDi)  (Eij) 

where  R^  j  is  the  reconstructed  value  for  the  i^  point  along 
the  profile  of  the  jth  eigenvector.  is  the  mean  profile 
elevation  at  the  ith  point  along  the  profile.  ±.C  is 
constant  which  is  representative  of  the  case  loadings  on  a 
particular  vector.  The  sign  is  dependent  upon  the  cartesian 
quadrant  being  reconstructed.  is  the  standard  deviation 

of  the  profile  elevation  at  the  ith  point  along  the  profile, 
and  E^j  is  the  loading  on  the  jth  eigenvector  for  the  ith 
point  along  the  profile  (Fig.  14  and  16) .  These  reconstruc¬ 
ted  profiles  for  each  quadrant  demonstrate  the  meaning  of 
various  points  in  eigenvector  space. 

The  passage  of  Hurricane  Dennis  (between  the  August  4 
and  August  24  surveys)  caused  substantial  changes  in  only 
two  FRF  profile  lines  (numbers  4  and  6).  Both  of  these 
profile  lines  were  located  on  the  north  side  of  the  pier. 
Profile  line  number  4  experienced  large  positive  movements 
while  line  number  6  showed  large  negative  movements.  Pro¬ 
file  line  numbers  3  and  6  were  the  only  profile  lines  to 
show  a  negative  movement  along  E2  during  the  period.  The 
remaining  profile  lines  showed  positive  movement  along  the 
E2  axis;  however,  most  of  these  were  small  increments. 

The  September  19  to  November  3  survey  period  showed 
smaller  changes  on  all  profiles  except  profile  line  number 


5.  These  changes  are  denoted  by  the  heavy  solid  lines  in 
Fig.  9  and  10.  All  profile  lines  except  number  11  showed 
large  negative  tendencies  in  the  E2  direction  with  very 
little  variance  in  El  values.  This  indicates  that  bar 
formation  occurs  prior  to  changes  in  the  slope. 

In  E1-E2  space  the  greatest  change  in  the  morphology 
occurred  between  the  November  3,  1981.  and  November  16. 
1981.  survey  dates;  and  all  profile  lines,  except  numbers  2 
and  8.  showed  these  changes.  This  change  is  shown  by  the 
heavy  dashed  black  lines  in  Fig.  9  and  10.  In  all  cases, 
the  weightings  on  E2  became  negative,  indicating  change  from 
quadrants  I  and  IV  to  quadrants  II  and  III.  A  negative 
weighting  on  E2  shows  a  seaward  displacement  of  sediment  and 
the  formation  of  a  double-bar  storm  profile.  During  the 
period,  loadings  on  El  generally  became  more  positive, 
indicating  offshore  displacement  of  sand.  Profile  line 
numbers  1-4  and  9-11  ended  in  quadrant  II  after  this  period 
with  a  well  developed  storm  bar  360  m  from  the  baseline  and 
a  smaller  inner  bar  140  m  from  the  baseline.  The  remaining 
profile  lines  (numbers  5-8).  which  are  located  in  close 
proximity  to  the  pier,  showed  a  well  developed  inner  bar  and 
a  smaller  storm  bar  located  in  the  same  locations  as  on  the 
other  profiles.  The  only  difference  in  the  profiles  was  the 
relative  size  of  the  two-bar  forms. 

The  plots  in  E3-E4  space  also  showed  large-scale 
variations  in  the  FRF  bathymetry  during  the  October-November 


period.  The  maximum  morphological  changes  on  all  profile 


lines  (except  numbers  4-6)  occurred  between  the  September  19 


and  November  3  surveys.  These  changes  are  shown  in  Fig. 


15a-15k  by  the  heavy  black  lines.  All  movements  were  into 


quadrant  II  where  a  storm  bar  is  centered  300  m  from  the 


baseline  and  a  small  inner  bar  is  welded  to  the  beach  face 


(Fig.  16).  These  changes  appear  to  designate  the  initial 


stages  of  the  storm  profile  development  and  the  offshore 


displacement  of  sediment  along  the  profile. 


The  largest  changes  observed  between  the  November  3 


and  the  November  16  surveys  occurred  in  profile  line  numbers 


4-7.  These  changes,  characterized  by  movements  into  quad¬ 


rant  If  indicate  a  net  loss  of  sediment  at  all  profile 


points  and  are  shown  by  the  heavy  dashed  lines  in  Fig. 


15a-15k.  Values  for  E4  have  large  positive  increments.  For 


this  same  period  Miller  et  al.  (1983)  described  large  gross 


changes  (the  sum  of  the  absolute  values  of  all  volumetric 


changes  along  a  profile)  away  from  the  pier,  while  net 


changes  (the  sum  of  all  positive  and  negative  elevational 


changes)  were  small.  In  the  vicinity  of  the  pier,  net 


changes  were  high.  They  noted  a  net  loss  of  sediment  from 


the  north  side  of  the  pier  to  the  south  side. 


Analysis  of  Shore-Parallel  Transects 


The  first  eigenvector  of  principal  components  analysis 


of  the  data,  organized  as  a  series  of  shore-parallel 
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transects/  explained  98%  of  the  total  variance.  The  first 
vector  indicates  either  more  or  less  sand  than  the  mean  at 
all  locations  along  the  shore-parallel  transect.  Because 
the  transects  are  shore  parallel,  variance  associated  with 
tars  and  troughs  is  not  present.  Only  relative  elevation 
between  points  is  shown.  None  of  the  higher  order  vectors 
met  the  Overland  and  Preisendorfer  (1982)  statistical  test/ 
nor  did  they  make  geophysical  sense.  This  analysis  was  also 
performed  on  the  data  separating  the  north  and  south  sides 
of  the  pier.  The  mean  loading  for  all  survey  dates  was 
calculated  for  the  jth  eigenvector  at  the  ith  point  along 
the  profile.  Each  of  the  44  mean  values  (E^)  was  then 
compared  to  the  actual  vector  loading  and  the  sign  of 
the  difference  was  plotted  in  Fig.  17.  providing  a  prognos¬ 
tic  overview  of  the  time  series  of  sediment  displacements. 
At  each  distance  from  the  shore  the  sum  of  E^t  over  all 
surveys  is  zero.  Thus,  the  chart  shows  the  variations  in 
the  sediment  budget  over  the  study  period. 

In  terms  of  a  gross  sand  budget  there  was  no  net 
aggregate  difference  for  one  side  of  the  pier  relative  to 
the  other.  More  sand  than  the  average  amount  was  in  the 
subaerial  region  during  surveys  1-6.  Sand  in  this  region 
shifted  to  less  than  the  average  amount  during  surveys  7-10. 
The  region  immediately  seaward  of  the  subaerial  region.  100 
to  175  m  from  the  baseline,  had  the  opposite  response.  The 
early  surveys  showed  less  sand  than  the  mean,  while  the 


■P  Oj  <1) 


B)  Hrl 


+*  *H  'O 


m  a)  cn 


1  ,  -C 

J-l  «H  -P 


al  0  3 


(0  <0  b  rH 

H'H  (3  ® 


[I] 


later  surveys  showed  more  than  the  mean.  We  concluded  that 
the  sand  was  shifted  from  the  subaerial  region  to  this 
region. 

Between  175  m  and  300  m  from  the  FRF  baseline  there 
was  a  relative  sand  surplus  in  the  early  surveys  and  a 
deficit  in  the  later  surveys.  The  apparent  recipient  of 
this  sand  in  the  later  surveys  was  the  zone  between  300  and 
600  m.  In  the  seaward-most  area,  600  to  800  m,  there  was, 
again,  a  sand  surplus  in  the  early  part  of  the  record  and  a 
deficit  in  later  periods.  While  the  profiles  do  not  extend 
seaward  of  800  m  from  the  FRF  baseline,  it  is  reasonable  to 
assume  that  the  "missing"  sand  in  the  later  surveys  was 
transported  to  some  point  beyond  800  m  from  the  FRF 
baseline.  The  resolution  of  the  data  in  this  region  makes 
the  whereabouts  of  the  "missing  sand"  somewhat  questionable, 
however. 

Most  of  the  observed  changes  occurred  between  the  5th 
and  the  7**1  surveys,  apparently  as  a  result  of  the  high 
pressure  "storms"  of  October  and  the  cyclones  of  November. 
It  is  interesting  that  the  exchanges  across  the  shoreline 
(100  m) ,  across  the  300-meter  distance,  and  across  some 
distance  seaward  of  800  m  occurred  synchronously.  The  time 
of  change  does  not  appear  to  be  a  function  of  water  depth. 
This  simultaneity  suggests  shore-normal  resonant  motions 
such  as  edge  waves.  If  this  is  the  case  then  this 
shore-parallel  form  of  eigenvector  analysis  may  prove  to 
have  been  especially  useful. 
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DISCUSSION  AND  CONCLUSIONS 


The  summer  to  winter  transition  in  shorezone  morphol¬ 
ogy  is  clearly  evident  in  the  principal  components  of  the 
correlation  matrix  of  the  FRF  shorezone  profiles.  During 
this  period  the  shorezone  became  more  dissipative  as  the 
profile  became  flatter  and  a  double-bar  configuration  re¬ 
placed  the  single  bar  of  summer.  These  changes  occurred  si¬ 
multaneously  over  all  water  depths.  In  terms  of  the  bar- 
and-trough  configuration  there  is  apparently  an  equilibrium 
form  "appropriate"  to  each  storm  that  occurs.  The  small 
differences  in  the  wave  characteristics  between  individual 
storms  were  insufficient  to  account  for  the  variations  in 
bar/trough  equilibrium  or  the  observed  changes  in  profile 
steepness. 

Two  hurricanes,  two  anticyclones,  and  two  cyclones 
occurred  during  the  period  of  study.  Each  of  these  storms 
generated  waves  of  similar  heights  (2  to  2.5  m) .  The  two 
hurricanes  produced  waves  that  resulted  in  nearly  identical 
bar/trough  equilibrium  forms.  Neither  hurricane  caused 
changes  in  profile  steepness.  The  two  anticyclones  of 
October  caused  the  profiles  to  flatten  and,  for  the  equilib¬ 
rium  bar/trough  configuration,  to  shift  markedly  toward  a 
two-bar  form.  The  two  cyclones  of  November  produced  waves 
that  caused  further  flattening  of  the  profile  and  accentua¬ 
tion  of  the  two-bar  configuration. 


The  waves  produced  during  the  storms  caused  adjust¬ 
ments  in  the  bar/trough  equilibrium  form;  whereas,  periods 
without  storms  showed  no  evidence  of  changes  in  the 
bar/trough  profile.  In  these  analyses  the  direction  of 
change  in  the  bar/trough  equilibrium  point  was  not  shown  to 
be  directly  related  to  either  wave  height,  wave  period,  wave 
direction,  wave  steepness,  breaker  type  or  tidal  levels, 
since  all  parameters  varied  very  little  between  storm 
periods.  The  highest  tide  during  this  study  occurred  in 
November  1981,  after  the  equilibrium  point  had  already  begun 
shifting  toward  the  two-bar  configuration.  It  was,  there¬ 
fore,  not  the  initiating  factor  in  the  profile  rotation. 
This  increase  in  the  mean  water  level  obviously  contributed 
to  the  overall  beach  changes,  however. 

The  only  attribute  of  the  stormy  periods  that  is 
possibly  associated  with  large  changes  in  profile  slope  and 
bar/trough  pattern  is  the  duration  of  unidirectional  long¬ 
shore  currents.  In  all  but  one  storm  event  (Tropical  Storm 
Emily),  southbound  longshore  currents  were  documented  for 
periods  of  eight  to  ten  days. 

Miller  et  al.  (1983)  have  documented  similar  situa¬ 
tions  at  the  FRF.  They  noted  that  the  pier  can  cause  very 
pronounced  effects  on  the  FRF  bathymetry.  The  shoreline 
updrift  of  the  pier  erodes,  while  it  builds  up  downdrift  of 
the  pier  during  periods  of  unidirectional  longshore  trans¬ 
port. 


Pre-existing  morphology  is  an  important  factor  in  con¬ 
trolling  the  subsequent  profile  changes  and  supports  the 
concept  of  sequential  beach  stage  models  proposed  by  Sonu 
and  Van  Beek  (1971)  and  Wright  et  al.  (1979).  At  the  Duck 
site  the  magnitude  of  change  toward  the  equilibrium 
bar/trough  configuration  is  a  function  of  the  bar/trough 
configuration  before  the  storm  occurred.  It  was  also  found 
that  the  magnitude  of  the  profile  slope  change  during  a 
storm  is  a  function  of  the  pre-storm  profile  slope. 

A  single  shoreline  crescentic  wave  is  usually  present 
at  the  Duck  site  with  a  wave  length  equal  to  the  reach 
surveyed  by  the  CRAB.  Crescentic  forms  in  this  size  class 
have  been  previously  reported  along  the  Outer  Banks  by 
Vincent  (1973).  The  amplitude  of  profile  elevations  along 
the  crescentic  form  is  greatest  seaward  of  300  ro  from  the 
FRF  baseline.  Because  the  crescentic  form  involves  the 
profile  rotation  eigenfunction,  it  seems  apparent  that  these 
crescentic  forms  are  directly  associated  with  onshore- 
offshore  sediment  exchanges.  During  the  October  and 
November  storms  sediment  was  transported  offshore. 

The  offshore  sediment  exchanges  of  October  and  Novem¬ 
ber  occurred  simultaneously  across  three  zones.  Sediment 
was  moved  from  the  beach  to  the  nearshore  where  a  bar  was 
formed.  Mid-profile  sediments  from  the  region  of  the  summer 
bar  were  moved  to  the  winter  bar  position.  Sediments 
seaward  of  800  m  were  moved  farther  seaward.  The  spacing 


between  these  three  zones  is  250  in.  The  November  changes 
in  the  subaerial,  mid-profile,  and  deep  water  occurred  at 
the  same  time  even  though  the  environmental  variables  were 
not  significantly  different  from  other  storm  events  where  no 
major  morphologic  changes  were  observed. 
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0112  2  ENVIRONMENTAL  RESEARCH  PRDCTN  FCLTY 

0112  3  NAVAL  POSTGRADUATE  SCHOOL _ 

0112  4  MONTEREY,  CALIFORNIA  93940 


0201  l  COMMANDING  GENERAL 

0201  2  MARINE  CORPS  DEVELOPMENT  AND 

0201  3  EDUCATIONAL  COMMAND 

0201  4  QUANTICO,  VIRGINIA  22134 


0202  1  OR.  A.  L.  SLAFKOSKY 

0202  2  SCIENTIFIC  ADVISOR 

0202  3  COMMANDANT  OF  THE  MARINE  CORPS _ J 

0202  4  CODE  MC-RD-1  . . . 

0202  5  WASHINGTON,  D.  C.  20380 _ ' 


1C01  1  DEFENSE  INTELLIGENCE  AGENCY 

1001  2  CENTRAL  REFERENCE  DIVISION 

1001  3  CODE  ROS-3 _ 

1001  4  WASHINGTON,  D.  C.  20301 


1002  1  DIRECTOR 

1002  2  DEFENSE  MAPPING  TOPOGRAPHIC  CENTER 

1002  3  ATTN=  CODE  50200  ; 

1C02  4  WASHINGTON,  D.  C.  20315  _ ! 

1003  1  COMMANDING  OFFICER 

1003  2  U.S.  ARMY  ENGINEERING 

1003  3  TOPOGRAPHIC  LABORATORY _ 

1003  4  ATTN=  ETL-ST 

1003  5  FORT  BELVOIR,  VIRGINIA  22060 


1005  1  CHIEF,  WAVE  DYNAMICS  DIVISION 

1005  2  USAE-WES 

1005  3  P.  C.  BOX  631 _ 

1005  4  VICKSBURG,  MISSISSIPPI  39180 


2C03  1  NATIONAL  OCEANOGRAPHIC  DATA 

2003  2  CENTER  ID764 _ 

2003  3  ENVIRONMENTAL  DATA  SERVICES 

2003  4  NOAA _ 

2003  5  WASHINGTON,  0.  C.  20235 


2005  1  CENTRAL  INTELLIGENCE  AGENCY 

2005  2  ATTENTION  OCR/DD-PUBL  IC ATI ONS 

2005  3  WASHINGTON,  D.  C.  20505 


SYSTEMATIC  GEGGRAPHY  LISTING 


i'vMH.- 


rjutumaM 


2007  2  ADVANCED  TECHNOLOGY  DIVISION 


uitw*na<H« 


(DELETED) 


2007  4  NAVAL  C3SERVAT0RY 


5001  4  HARDTHGEHE 

0 


5002  2  RUE/FO 


500  3 


5C03  2 

500  3  3 


5003  - 


5302  1 

5302 


5302  3 

5302 _ 4 


DEFENCE  RESEARCH  ESTABLISHMENT 
CSTER9R0GADES  KA 


DK-2100  KCBENHAVN  0,  DENMARK 


IR.  M.  Vi.  VAN  BATENBERG 
PHYSISCH  LABORATCJRIUM  TNO 


UUDE  WAALSDORPER  REG  63 ,  DEN  HAAG 
N  ETHERLAND 


8002  2 
8002  3 


LOUISIANA  STATE  UNIVERSITY 
BATGN  ROUGE.  LOUISIANA  70 


8003  4  SIERRA  MADRE,  CA  91024 


9C01  1  DR.  LESTER  A.  GERHARDT 

RENNSSELAFR  POLY  TECH 


9001  3  TROY.  NEW  YORK  1218  1 


9002  1  MR.  FRED  THOMSON 

M' 


U'RVattlll 


9002  3  P.G.  BOX  618 


CCT.  26,  1977 


JOB  460-101 


900  4 

1 

SYSTEMATIC  GEOGRAPHY  LISTING 

DR.  THOMAS  K.  PEUCKER 

9004 

2 

SIMCN  FRASER  UNIVERSITY 

9004 

3 

DEPARTMENT  OF  GEOGRAPHY 

9C04 

4 

BURNABY  2 ,  B.  C.,  CANADA 

9005 

1 

DR.  ROBERT  DCLAN 

9005 

2 

DEPARTMENT  OF  ENVIRONMENTAL  SCIENCES 

9005 

3 

UNIVERSITY  OF  VIRGINIA 

9005 

4 

CHARLOTTESVILLE,  VIRGINIA  22903 

